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Abstract 

We present the matching condition for the strong conpling con- 
tant at a heavy quark threshold to four loops in the modihed 
minimal subtraction scheme. Our results lead to further decrease 
of the theoretical uncertainty of the evolution of the strong cou¬ 
pling constant through heavy quark thresholds. Using a low en¬ 
ergy theorem we furthermore derive the effective coupling of the 
Higgs boson to gluons (induced by a virtual heavy quark) in four- 
and (partially) through hve-loop approximation. 
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1 Introduction 

The masses of the known quark species differ vastly in their magnitude. As a 
result of this, in many QCD applications the mass of a heavy quark h is much 
larger than the characteristic momentum scale ^/s which is intrinsic to the 
physical process. In such a situation there appear two interrelated problems 
when using an MS-like renormalization scheme. 

First, given the two large but in general quite different mass scales, ^/s 
and m, two different types of potentially dangerously large logarithms may 
arise. The standard trick of a proper choice of the renormalization scale fi 
is no longer effective; one can not set one parameter fi equal to two very 
different mass scales simultaneously. 

Second, according to the Appelquist-Carazzone theorem [1] the effects due 
to heavy particles should eventually “decouple” from the low-energy physics^. 
However, a peculiarity of mass-independent renormalization schemes is that 
the decoupling theorem does not hold in its naive form for theories renor¬ 
malized in this framework. The effective QCD action to appear will not be 
canonically normalized. Even worse, potentially large mass logarithms in gen¬ 
eral appear, when one calculates a physical observable. 

The well-known way to bring the large mass logarithms under control is 
to construct an effective field theory by “integrating out” the heavy quark 
field h [2-6]. By construction, the resulting effective QCD Lagrangian will 
not include the heavy quark field. In order to be specific, let us consider 
QCD with ni = Uf — 1 light quarks and one heavy quark h with mass m. 
The quark-gluon coupling constants in both theories, the full n/-flavor QCD 
and the effective nz-flavor one, and are related by the so-called 

matching condition of the form [2, 3] 

( 1 ) 

Here m = m{fi) is the (running) mass of the heavy quark; the decoupling 
function Q is equal to one at the leading (tree) level but receives nontrivial 

^It should be stressed that the statement is literally valid only if power-suppressed 
corrections of order (s/m^)" with n > 0 are neglected. 
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corrections in higher orders. The matching point /i in eq. (1) should be chosen 
in such a way to minimize the effects of logarithms of the heavy quark mass, 
e.g. /i = 0{m). 

An important phenomenological application of eq. (1) appears in the 
determination of as{Mz) at the Z-boson scale through evolution with the 
renormalization group equation, starting from the measured value as{mr) 
at the r-lepton scale. A careful analysis of the effects of four-loop running 
and three-loop matching in extracting af\Mz) with 5 active quark flavors 
from a^pirrir) with 3 active quark flavors has been recently performed in 
ref. [7]. We will demonstrate that the inclusion of the newly computed four- 
loop matching condition leads to further reduction of the theoretical error 
from the evolution. 

The second aspect, mentioned above, is of importance for Higgs-boson 
production in hadronic collisions. The dominant subprocess for the produc¬ 
tion of the Standard-Model (SM) Higgs boson at the CERN Large Hadron 
Collider (LHC) will be the one via gluon fusion. Therefore, an important in¬ 
gredient for the Higgs-boson search will be the effective coupling of the Higgs 
boson to gluons, usually called Ci. In ref. [8] a low-energy theorem, valid to 
all orders, was established, which relates the effective Higgs-boson-gluon cou¬ 
pling, induced by the virtual presence of a heavy quark, to the logarithmic 
derivative of C,g with respect to the heavy quark mass. In that paper the 
method was used to squeeze from the three-loop decoupling function C,g the 
analytical result for Ci in three- and even in four-loop approximation (the 
latter in a indirect way through a sophisticated use of the RG evolution equa¬ 
tions and the four-loop QCD /5-function). With our new full four-loop result 
for Qg we are able to confirm the result of ref. [8] in a completely independent 
way, without using the four-loop contribution to the QCD /5-function. In ad¬ 
dition a five-loop prediction for Ci can be obtained (modulo yet unknown 
contribution from the five-loop nj-dependent term in the /5-function). 

The outline of the paper is as follows. In Section 2 we recall the main 
formulae from [8] which reduce the evaluation of the decoupling function C,g 
to the calculation of vacuum integrals. Then we discuss shortly the technique 
used to compute these integrals. Section 3 describes our four-loop results for 
the decoupling function. In Section 4 we make use of the low energy theorem 
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to derive the effective coupling of the Higgs boson to gluons (due to the 
virtual presence of a heavy quark) through four and (partially) through hve 
loops. Section 5 deals with phenomenological applications of our result which 
lead to a reduction of the theoretical error due to evolution of the coupling 
constant through heavy quark thresholds. Summary and conclusions are 
presented in the hnal Section 6. 


2 Formalism 


2.1 Decoupling for the Gauge Coupling Constant 

Let us consider QCD with ni = nj — 1 massless quarks tjj = = l-n^} and 

one heavy quark h with mass m. The corresponding bare QCD Lagrangian 
reads 

C{go, mo, ^o; ^o, c°’“, ho) = -^(<^2’“)^ + #o Wo + ho{ipo - mo)ho 

+ terms with ghost fields and the gauge-fixing term, 

( 2 ) 

where the index ‘0’ marks bare quantities, and are the gluon and 
ghost held respectively. The relations between the bare and renormalized 
quantities read 


9s = = Zmm, f -1 = ^ 3 (^ - 1 ), 

y“=yi2<*„ = c“-‘‘ = \/Tc“, (3) 


where gs = y/^nag is the (renormalized) QCD gauge coupling, g, is the renor¬ 
malization scale, d = 4 — 2e is the dimensionality of space-time. The gauge 
parameter ^ is dehned through the gluon propagator in lowest order. 


i 

+ ie 


- 9 ^^ + e 




( 4 ) 


All renormalization constants appearing in (3) are known by now through 
order 0(ag) from [9-13]. 

Integrating out the heavy quark transforms the full QCD Lagrangian 
(2) into the one corresponding to the effective massless QCD with n'j = 
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Uf — 1 = ni quark flavors (plus additional higher dimension interaction terms 
suppressed by powers of the heavy mass and neglected in what follows). 
Denoting the effective helds and parameters by an extra prime, we write the 
effective Lagrangian as follows: 

= (5) 

Here the primed quantities are related to non-primed ones through 

j,"'=o”. ?“'-!= - 1 ). <= 

G“''“ = ( 6 ) 


where the primes mark the quantities of the effective nrflavor theory. 

As was hrst demonstrated in [8] the bare decoupling constants can all be 
expressed in terms of massive Feynman integrals without any external mo¬ 
menta (so-called massive tadpoles). The corresponding relation for (jj reads: 



where 


(7) 


i + r"t,(o,o). 

(8) 

i + n^"(o), 

(9) 

i + nf(o). 

(10) 


Here UgIp^) and are the gluon and ghost vacuum polarization func¬ 
tions, respectively. and are related to the gluon and ghost 

propagators through 


^ab 


fiu 


y [i + n^(p2); 


-|- terms proportional to p^p‘' 


^ab 

[1 + no(p2)] 


i I dxe^P-^ (TG'°’“^(a;)G'°’^"(0)), 
i I (Tc^’dx)d>’\0)) , 


( 11 ) 


( 12 ) 
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respectively. The vertex function V^Q^{p, k) is defined through the one-particle- 
irreducible (IPI) part of the amputated Gcc Green function as 




|_^ya6c k) + other color structuresj 


= e dxdy (Tc^’%x)^’\0)G^’^y{y) 


ipi 


(13) 


where p and k are the outgoing four-momenta of c and G, respectively, and 
are the structure constants of the QCD gauge group. 

Finally, in order to renormalize the bare decoupling constant (jj we com¬ 
bine eq. (7) with eqs. (3) and arrive at (q;^ = 5f^/(47r),= (5f()^/(47r)) 


a',(y) 



Cg*^s(h')- 


(14) 


Direct application of this equation is rather clumsy as the constant Z' on its 
right side depends on the renormalized effective coupling constant which we 
are looking for. Of course, within perturbation theory, one can always solve 
eq. (14) by iteration. A simpler way, which we have used, reduces to inverting 
hrst the series 


a® — Za(as, e) as, Za{as, e) — Z^ — 1 + ^ Za^i{€) a* (15) 

i>l 

to express the renormalized coupling constant in terms of the bare one: 

a, = Z°(a°,e)a°, Z°(a„ e) = 1 + ^ Zj.(e) (a°)\ (16) 

i>l 


where, up to four loops. 


vu — 
^0,1 ~ 

70 _ 

^ 0,4 — 


ctPi ^a.2 ~ ~^oi,2 1 " ^ -^ 0,3 — ~Za,3 + 5 Z^^l Za^2 

2 oi ^2 V I 1 4 ^4 (17) 


-Zn3 , Z„, o — ~Zn 2 2 Z^ 

— Za,4: + 6 Za,l Za^3 + 3 Z‘^ 2 ~ 21 Z^ -^ Za^2 + 14 Z^^^ 


5^1 


With the use of eq. (17) one now could conveniently transform all the primed 
(that is effective quantities ) in eq. (14) into the non-primed ones. After this 
the renormalization can be done directly. 


5 



2.2 Vacuum Integrals 

At the end of the day we have to evaluate a host of four-loop massive tadpoles 
entering the dehnitions of the bare decoupling constants. The evaluation of 
these massive tadpoles in three-loop approximation has been pioneered in 
ref. [14] and automated in ref. [15]. 

Similar to the three-loop case, the analytical evaluation of four-loop tadpole 
integrals is based on the traditional Integration-By-Parts (IBP) method. In 
contrast to the three-loop case the manual construction of algorithms to re¬ 
duce arbitrary diagrams to a small set of master integrals is replaced by 
Laporta’s algorithm [16,17]. In this context the IBP identities are generated 
with numerical values for the powers of the propagators and the irreducible 
scalar products. In the next step, the resulting system of linear equations is 
solved by expressing systematically complicated integrals in terms of simpler 
ones. The resulting solutions are then substituted into all the other equa¬ 
tions. 

This reduction has been implemented in an automated FORMS [18,19] based 
program in which partially ideas described in ref. [17,20, 21] have been imple¬ 
mented. The rational functions in the space-time dimension d, which arise in 
this procedure, are simplified with the program FERMAT [22]. The automated 
exploitation of all symmetries of the diagrams by reshuffling the powers of 
the propagators of a given topology in a unique way strongly reduces the 
number of equations which need to be solved. 

In general, the tadpole diagrams contributing to the decoupling constants 
contain both massive and massless lines. In contrast, the computation of the 
four-loop /5-function can be reduced to the evaluation of four-loop tadpoles 
composed of completely massive propagators. These special cases have been 
considered in [9,13,21]. 

All four-loop tadpole diagrams encountered during our calculations were 
expressed through the set of 13 master integrals shown in figures 1 and 2. 
The first four integrals of Fig. 1 have analytical solution in terms of Gamma 
functions for generic valuse of e. The first nine terms of the e-expansion of 
the hfth integral (T 52 ) can be obtained from results of Refs. [14,23] (for 
more details see Ref. [28]). As for the less simple master integrals pictured 
in Fig. 2 one hnds that for the case of (g all necessary ingredients are known 
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analytically ([24-29]) except for T 54 ^ 3 , Tq 2,2 and Tgi^, where we have denoted 


T ■ 

n 



d^-^-p 1 _ 

^2-e (jp, _|_ 1-^2 J ^^0 


E 

i>-4 


(18) 


These three integrals are known numerically from [28,29] and read 


^ 54,3 = 

-8445.8046390310298, 

(19) 

^62,2 = 

-4553.4004372195263, 

( 20 ) 

Tgip = 

1.808879546208. 

( 21 ) 


In fact, in refs. [29, 30] an analytical result for Tgi^o (in terms of Te 2,2 and 
^ 54 , 3 ) has been also derived. Thus, the results of the next section will contain 
only two numerical constants, viz. Tq 2,2 and 



Figure 1: Analytically known master integrals 


3 Decoupling at four loops: results 

We have generated the relevant diagrams with the help of the program 
QGRAF [31]. The total number of diagrams at four loops amount to 6070, 
765 and 9907 for and Ci respectively. All calculations were performed in 
a general covariant gauge with the gluon propagator as given in eq. (4). How¬ 
ever, since extra squared propagators lead to signihcant calculational compli¬ 
cations, only linear terms in ^ were kept. The bare results are rather lengthy 
and will be made available (in computer-readable form) in http://www- 
ttp.physik. uni-karlsruhe. de/Progdata/ttp05 / ttp05-27. 

After constructing the bare decoupling constant from eq. (7) followed by 
the renormalization as described in section 2 we arrive at the following gauge 
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T54,(4,3,2) 


Ts2. (3,2,1) 


T8„(-2,1,1) 


Ts4, (2.1,1) 



T„,(1,0,0) 



T72, (1, 0, 0) 




T8„(0,-1,-1) T9„(1,0,-1) 


Figure 2: Master integrals where only a few terms of their e-expansion are 
known analytically. The solid (dashed) lines denote massive (massless) prop¬ 
agators. The three numbers in brackets {ni,n 2 ,n^) are decoded as follows: 
ni is the maximal power of the spurious pole in e which might appear in 
front of the integral, n 2 is the maximal power of the spurious pole in e which 
happens to enter into the decomposition of the bare deconpling constant (jj 
in terms of the master integrals, is the maximal analytically known power 
of the e-expansion of the same integral as determined in [28] and conhrmed 
in [29]. 


independent^ result for the decoupling function 

(gMS)2 ^ dMS,i, 


i>l 


where we use the notation 






TT 


and the coefficients dus,! read 


( 22 ) 


dus,! — — 7 : 


g -fim 1 


(23) 


^In fact, the gauge dependence on ^ disappears already for the bare decoupling constant 
as it should be. 







( 24 ) 


_11 11 1 2 

“Ms ,2 — ^ 36^^™’ 


<^MS,3 — 

+ 


564731 


124416 

ni 


82043 955 

Ca — 


27648 
2633 67 

+ 


576 

1 , 


53 ,2 

-I-7^ 

576 


216 


31104 576 


_ )2 

36 


(25) 


C^MS,4 — ^MS,4 + 


+ ni 


7391699 2529743 

746496 ~ 165888 

110341 110779 

373248 000/14 ’3 iim 


+ 2177 ,, 


1883 


3456 10368 

1483 127 


+ 


1 


+ n. 


6865 
.186624 


82944 
77 „o 


^(im 

1 


10368 5184 


_ I ~ p'i 

20736^^^ 324 


1296 


(26) 


'^MS,4 — 


134805853579559 254709337 


TT^ — 


151369 


43342154956800 783820800 30481920 

4330717 ^ 9869857 121 2057 

Cs H-®4 — -— 


TT^ - 


18233772727 , 151369 

j'J r A A oor\ ^3 


783820800 


544320 


207360 


36 


51840 


272160 

121 2, a 9869857, . 121 , . 

TT^ In^ 2 +-In^ 2 +-In® 2 + 

' c-\ -t A ' A c-\ f~\ ' 


7rMn2 - 


9869857 


TT^ In^ 2 


6531840 
82037 „ 151369 


+ 


ni 




2592 
4770941 

"2239488 

685 

~ 124416 
271883 


■T; 


54,a 


6531840 4320 30965760 

541549 4 3645913 115 685 

Ca + Cs + TTTTT ®4 


11612160 


■Tf 


62,2 


14929920 
7r2ln2 2 + 


995328 


576 


5184 


685 


124416 


-ln^2 


+ 


167 


4478976 5184 


Ca 


(27) 


In eqs. (23-27) = In ^ 2 (^) , rn{^) is the (running) heavy quark mass 

in the MS-scheme and /i represents the renormalization scale. Furthermore, 

Cn = C(^) is Riemann’s zeta function and = LR(l/2) = l/(2*i"'). 

For another convenient dehnition of the quark mass — the so-called scale- 

invariant mass, dehned through the relation Hh = m{fih) — eqs. (22-27) are 

2 

transformed to {i^h = In ^): 


(cfy = 1 +!!<(>■) cisi.i, 


(28) 


i>l 
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where 


<^si,i — — -p: ^^ih , 
0 


11 19 1 

dsi,2 “ 72 “ ^ ^ 


564731 82043 6793 

“ 124416 ~ 27648 ~ 17^ 

\ 2633 , 281 ^ ] 


_^73 

576 216 


31104 1728 


^ 2398621 2483663 14023 ^ 8371 3 1 4 

SI,4 - MS,4 - 74g49g - ig5333 ^3 ^^lh “ 

^90283 133819 983 2 107 3 ' 

^ [373248 ^ 82944 ^ 3456 ^ 

2 r 8545 „ 79 ,2 1 

^ [ 186624 ~ Qm 2 ^ ^ 


• practical applicatl 


Ions also the inverted formulae are needed: 
1 . ... , 


/^MS'|2 - 1 (o^MT d^ 

\^g } i>l 

77^ = 1 + H («s(/^))* 4 

v8g ) i>l 


where we use the notation 

aM = 

IT 

and the coefficients dj^g;, dgj j read 


d' —- i 

“MS,1 ~ g 1 


, _ _11 11 ^ ]_ p ‘2 

'^MS,2 ~ 72 24 ^ 36 
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'^MS,3 ~ 


564731 82043 , 2645 


124416 ^ 27648 
2633 67 


Cs + 


„ I . 

- — - ^iim 

[31104 576 ^36 ' 


1728 

— f 

36 


167 ,2 1 .3 

H- r A - 

576 216 


( 37 ) 


d' — 
“ms,4 ~ 


121 11093717 

^MS,4 


ni 


+ nj 


1728 

141937 


3022001 , , 

- ^ nm 4 “ - C 3 ^ 

746496 ^ ICKOOO , 


373248 
6865 


110779 
82944 


+ 


165888 
277 


7^m 


f2 271 3 

10 ^ Vm + ^ 

1 -• 


6865 , 77 _ J_ 

186624 20736 324 


1837 ^2 2909 ^3 1 ^4 

J- -\ - P -\ - 

1152 10368 1296 

L 1 


(38) 




(39) 


d' — 
“SI,3 “ 


d' 

“SI,4 


si ,2 ^2 ^ 24 ^ 36 

564731 82043 2191 511 2 

124416 ^ 27648 ^ ^7^ ^ ^ ^ 

r 2633 281 

^ [31104 ~ 17^ ’ 


(40) 


_ 

216 


(41) 


121 ^ 1531493 2975921 33887 2 14149 3 1 4 

1728 ~ ~ 746496 165888 ^ 3456 ^ 10368 1^ 

158687 133819 515 2 107 3 ' 

'373248 ~ 82944 ~ 1152 ~ 

8545 


ni 


n, 


186624 ^ 6912 

'Numerically, eqs. ( 22 ) and (28) read 

(cy- 


(42) 


fl-- 


, _ ^-y - ^-y- 

= 1 + 0.152778 a2(/i;,)+ a3(/r;,) (0.972057-0.0846515 nd 

-fj-h 

+ at(/ift) (5.17035 - 1.00993 n«- 0.0219784 (4 
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1 




0.152778af(/ife) + a'/(/i/,) (-0.972057 + 0.0846515 m) 



1 - 


+ (-5.10032 + 1.00993 m + 0.0219784 . (44) 


Using the three loop relation between the pole and MS masses [32-34] 
one could easily express the decoupling relations in terms of the pole mass 
of the heavy quark. As the resulting expressions are rather lengthy they are 
relegated to the Appendix. 


4 Coupling of the Higgs boson to gluons 

Within the Standard Model (SM) the scalar Higgs boson is responsible for 
the mechanism of the mass generation. Its future (non-)discovery will be of 
primary importance for all the particle physics. The SM Higgs boson mass 
is constrained from below, Mh > 114GeV, by experiments at LEP and SLC 
[35, 36]. Indirect constraints coming from precision electroweak measurements 
[37] lead to an upper limit of about 200 GeV. 

With the SM Higgs boson mass within this range, its coupling to a pair 
of gluons is mediated by virtual quarks [38] and plays a crucial role in Higgs 
phenomenology. Indeed, with Yukawa couplings of the Higgs boson to quarks 
being proportional to the respective quark masses, the ggH coupling of the 
SM is essentially generated by the top quark only. The ggH coupling strength 
becomes independent of the top quark mass Mt in the limit Mh ^Mt. 

In general, the theoretical description of such interactions is very com¬ 
plicated because there are two different mass scales involved, Mh and Mt. 
However, in the limit Mh ‘^Mt, the situation may be greatly simplihed by 
integrating out the top quark, i.e. by constructing a heavy-top-quark effective 
Lagrangian [39,40]. 

This Lagrangian is a linear combination of certain dimension-four op¬ 
erators acting in QGD with Eve quark flavors, while all Mt dependence is 
contained in the coefficient functions. The hnal renormalized version of £eff 
is {Gp is Fermi’s constant) 

/leff = [0(]. (45) 


12 



Here, [O'J is the renormalized counterpart of the bare operator O'^ = 
where Ga^j.u is the color field strength, the superscript 0 denotes bare helds, 
and primed objects refer to the hve-flavor effective theory. Ci is the corre¬ 
sponding renormalized coefficient function, which carries all Mt dependence. 

In ref. [8] a low-energy theorem was established, valid to all orders, which 
relates the effective coupling of the Higgs boson to gluons, induced by a 
presence of heavy quark h, to the logarithmic derivative of C,g w.r.t. m. The 
theorem states that 

Cl = -]-m^ (46) 

2 om^ ^ 

Another, equivalent, but more convenient form of (46) was also derived in [8] 
by exploiting evolution equations in full and effective theories. It reads 


C'i = 


TT 


2 [1 - 2'ym{as)'. 


P'ia's) - P{as) - P{as)as-^ In Cg 


(47) 


where is the quark mass anomalous dimension, (g = Cg (h) as, m) and the 
/5-function is dehned as follows: 


yAot'’ = y”'> (oy) 


d /i 


OO 


N=1 


n{nf) N+l 


(48) 


a, = 


= al"'). 


asPias) = P^^^\as), a{P'{a{) = P^^^\a{). 


An important feature of the low energy theorem (47) is the fact that the 
decoupling function (^g appears there multiplied by at least one power of a*. 
It means that in order to compute Ci , say, at four loops one should know 
the QCD /5-function to four loops ( only Uf dependent part) but the quark 
anomalous dimension and the decoupling function (g only to three loops. 
It is this fact which allowed to hnd Ci at four loops in [8] long before the 
four-loop result for the decoupling function would be available. 

Now, armed with the newly computed four-loop term in (g, we could 
easily check that the old result for Ci of [8] by a direct use of eq. (46). We 
have done this simple exercise and found full agreement. In fact, one could 
even use eq. (47) in order to construct Ci at hve loops in terms of the known 
four-loop QCD decoupling function, the quark anomalous dimension y^ and 
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the /9-function and the only yet unknown parameter: the n/-dependent piece 
of the hve-loop contribution to the /9-function. The result reads 


C'i = - 

+ 


12 




2821 


+ ni[- 


67\ 


4004351 1305893,,^, 

+ C(3)+n; 


62208 


13824 


96/ 

/115607 110779 

V 62208 ~ 13824 


C(3)^ + 


6865 \ 
31104/ 


13546105 91 

41472 Y^ms,4 


31 


432 


■ TT 


+ ni 


+ ni 


/28598581 
V 497664 
3503 


3843215 


29 

3 MS,4 - Y2^ ^ 110592 


853091 475 

6912 ^ 

575 


Cs —^ Cs 

36 


1 4 3 ^ \ 31 

62208 ^216^ 4^3J + r^/| 


83 

7776 


— Cs 1 + 6 A//4 


12 


(h/; 


1 + 2.7500 as{^^h) + (9.7951 - 0.6979 ni) 

+ (49.1827- 7.7743 n/ - 0.2207nf) 

+ (-662.507 + 137.601 m - 2.53666 n^i - 0.077522 nf + 6 A/ 94 ) at(/ih) 


(49) 


where 

A/94 = /9f - /9f 

To save space we have written eq. (49) for the value of /i = fXh- Unfortu¬ 
nately, simple estimations show that the /9-function dependent contribution 
to eq. (49) could be numerically important. 


5 Application: as{Mz) from as{M^) 

A central feature of QCD is the possibility to describe measurements per¬ 
formed at very different energy scales with the help of only one coupling 
constant. The customary choice is the (MS) coupling constant as{Mz) which 
can be measured precisely in Z-boson decays. (For a review see, e.g. [41].) 

On the other hand the dependence of the r-decay rate on the strong 
coupling as has been used for a determination of at lower energies, with the 
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most recent results of 0.340 ±0.005exp± 0.014th and 0.348 ±0.009exp±0.019th 
by the ALEPH [42] and OPAL [43] collaborations. 

After evolution up to higher energies (taking into account the threshold 
effects due to the c- and 6-quarks) these results agree remarkably well with 
determinations based on the hadronic Z decay rate, which provides a genuine 
test of asymptotic freedom of QCD. 

Very recently an analysis of the evolution has been performed in [7]. They 
start from an updated determination of asirrir) 

as{mr) = 0.345 ± 0.004exp ± 0.008th (50) 

based on most recent experimental results of the ALEPH collaboration [42]. 
Their result for the evolution of asirrir) given in eq. (50), based on the use 
of the four-loop running and and three-loop quark-flavor matching reads 

asiMz) = 0.1215 ±0.0004exp± 0.0010th ± 0.0005evoi , 

= 0.1215 ± 0.0012 , (51) 

where the hrst two errors originate from the asirrir) determination given 
in Eq. (50), and the last error stands for the ambiguities in the evolution 
due to uncertainties in the matching scales of the quark thresholds. That 
evolution error received contributions from the uncertainties in the c-quark 
mass (0.00020, /ic = 1-31 ± 0.1 GeV) and the 6-quark mass (0.00005, fib = 
4.13 ± 0.1 GeV), the matching scale (0.00023, /i varied between 0.7/Xg and 
3.0 Hq), the three-loop truncation in the matching expansion (0.00026) and 
the four-loop truncation in the RGE equation (0.00031). (Eor the last two 
errors the size of the shift due the highest known perturbative term was 
treated as systematic uncertainty.) The errors had been added in quadrature. 

We have repeated this analysis including the newly computed four-loop 
approximation for the matching."^ As a result the value of asiMz) has been 
marginally increased (by 0.0001) and both errors from the choice of the 
matching scale and from the four-loop truncation in the matching equation 
have been halved. The updated version of (51) now reads 

asiMz) = 0.1216 ±0.0004exp± 0.0010th ± 0.0004evoi , 

= 0.1216 ± 0.0012 . (52) 

^We have used the package RunDec [44]. 
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6 Summary and Conclusions 


We have computed the decoupling relation for the QCD quark gluon cou¬ 
pling constant in four-loop approximation. The new contribution leads to a 
decrease of the matching related uncertainties in the process of the evolution 
of the as{mr) to as{Mz) by a factor of two. 

As a by-product we have directly conhrmed the long available result for 
Cl, the effective coupling of the Higgs boson to gluons at four loops, and (par¬ 
tially) extended it to hve loops. It remains only to hnd the QCD /5-fnnction at 
five loops to get the full result for Ci. In the light of recent signihcant progress 
in the mnltiloop technology [45,46] the completely analytical evaluation of 
the former seems to be possible in not too far distant future. 

We would like to mention that the result for the decoupling function at 
four loops as well as for the Higgs boson to gluons effective coupling at five 
loops have been obtained by a completely independent calculation in [30]. 
The results are in full mutual agreement. 
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A (g for the heavy quark mass in the on-shell 
scheme 


The relevant formulas for the case of the heavy quark mass renormalized in 
the on-shell scheme (denoted as M) are: 

(C°v = l + !:“•>) ‘ios.i (53) 

i>l 

and ^ 

7Tos\2 ^ ^ + X! (^'oS,v (54) 

\^g ) i>l 

2 

where if,M = In and 

dos,i ^ ~ Q ) (55) 


dos ,2 — 


24 


19 


^liM + 



(56) 


dos,3 


58933 1 _2 

~ 124416 ~ 9^ 
_J_^3 
576 216 


80507 , 1 2 , 

-Cs-TT^ In 2 

27648 27 

r 2479 1 2 

^ [31104 


8521 


1728 

409 1 

+ ^728 , 


, (57) 


dos,A 


^OS,4 

7693 


+ ni 


+ nf 


1152 

1110443 


19696909 29 2 . 

746496 M ^ 

10368 1296 


2439119 29 

165888 


373248 ^ 

1728 

1679 

186624 


41 2 « 132283 


TT^ In 2 4 m 


-'-2 1 /, 6661 2 

g^Cs^Af + sT’" >"2 4 m + Y03^4m 


-TT^ iaM 

162 ^ 


493 


n2 

20736 


(58) 


17 



^OS,4 — 


2180918653146841 


43342154956800 116640 


697121 2 231357337 . 

TT — _ TT — 


783820800 


151369 6 18646246327 

TT — 


30481920 


783820800 


Cs 


ni 


+ n, 


1439 2 . 151369 2 3698717 10609057 121 

^ TT A A nr\r\ ^3 


1296 

1027 

121 


TT^ In 2 


4320 
1773073 


In® 2 + 


544320 

2057 

' 51840 
82037 


207360 


TT^ In 2 


30965760 


T. 


54,3 


9278497 

' 6531840 
151369 

11612160 


272160 


04 — 


36 


05 


121 2 , 3 10609057, 4 

■ TT^ In® 2 + _ In^ 2 


2592 


Tf 


62,2 


557 


- 1 -TT^- 

746496 972 14929920 


697709 4 4756441 , 115 , 

Tt" + C3 + T— C5 


173 


04 


5184 
140825 

' 1492992 


11 


7rMn2 - 


13 2 19 

~ 1728 


1709 

124416 

C 3 


995328 
7r2ln2 2 + 


576 


173 

124416' 


-ln ^2 


6531840 


(59) 


'^OS,l —Q , 


(60) 


“os ,2 “ 24 24 36 


(61) 


'^OS,3 “ 


58933 


d' — 
“os,4 “ 


124416 9 


1 2 

+ 75 75 ^ + 


80507 ^ 1 21 ^ 8941 ^ 

27648 ^ ^ ^ ^ 


511 




576 216 

r 2479 1 2 

~ 7777777 ~ TT 


31104 54 


409 

17^ 
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(62) 


ni 


n, 


21084715 35 2 . 2922161 

^ 746496 + 54 ^ ^ 165888 

35 2 , 0 . 47039 ,2 14149 1 ,4 

162 ^ 3456 10368 1296 

1140191, 47 2 . 132283 ^, 1 2 . 

373248 ^ 324 ^ 82944 ^ 81 


ij 

'fiM 


107 

- ^728 
■ 1679 1 2 , 

+Y^7r + 


9115 

10368 


493 




(63) 


18 



((Osy ^ ^ _ 0.291667a2(M)+ a^(M) (-5.32389 + 0.262471 nz) 

^ )i=M 

+ a^(M) (-85.875 + 9.69229 m - 0.239542 , (64) 

^ 1 + 0.291667<2(M)+ a(3(M) (5.32389- 0.262471 nz) 

((■ 03)2 s \ / s \ \ / 

+ a(^(M) (86.1302 - 9.69229 m + 0.239542 . (65) 
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